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Abstract. We study a network of coupled logistic maps whose interactions occur with a certain distribution 
of delay times. The local dynamics is chaotic in the absence of coupling and thus the network is a paradigm 
of a complex system. There are two regimes of synchronization, depending on the distribution of delays: 
when the delays are sufficiently heterogeneous the network synchronizes on a steady-state (that is unstable 
for the uncoupled maps); when the delays are homogeneous, it synchronizes in a time-dependent state 
(that is either periodic or chaotic). Using two global indicators we quantify the synchronizability on the 
two regimes, focusing on the roles of the network connectivity and the topology. The connectivity is 
measured in terms of the average number of links per node, and we consider various topologies (scale-free, 
small- world, star, and nearest- neighbor with and without a central hub). With weak connectivity and weak 
coupling strength, the network displays an irregular oscillatory dynamics that is largely independent of 
the topology and of the delay distribution. With heterogeneous delays, we find a threshold connectivity 
level below which the network does not synchronize, regardless of the network size. This minimum average 
number of neighbors seems to be independent of the delay distribution. We also analyze the effect of self- 
feedback loops and find that they have an impact on the synchronizability of small networks with large 
coupling strengths. The influence of feedback, enhancing or degrading synchronization, depends on the 
topology and on the distribution of delays. 

PACS. 5.45.Xt, 05.45.-a, 05.45. Ra, 05.45. Pq 
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> . 1 Introduction 

, Complex networks arise due to self-organization plienoni- 
■ - - I ena in many real systems, such as food webs, the Inter- 
net, social networks, genes, cells and neurons [T]. A lot 
of research has been devoted to understanding the collec- 
tive behavior emerging in complex networks [5] , given the 
individual dynamics of the nodes and the coupling archi- 
tecture. Ecological webs, for example, describe species by 
means of nodes connected by links, representing the inter- 
actions. The interactions can be either direct or indirect 
through intermediate species, and can be of antagonistic 
type, such as predation, parasitism, etc., or mutually ben- 
eficial, such as those involving the pollination of flowers 
by insects |3] . A study of the network architecture decom- 
pose food webs in spanning trees and loop- forming links, 
revealing common principles underlying the organization 
of different ecosystems U . Cells use metabolic networks of 
interacting molecular components in processes that gener- 
ate mass, energy, information transfer and cell-fate speci- 
fication [S]. Adaptation and robustness have been shown 
to be consequences of the network's connectivity and do 
not require the 'fine-tuning' of parameters [5]. In a cat's 



brain, functional connectivity has been studied within the 
framework of complex networks [7]- In the human brain, 
magnetic resonance imaging has been used to extract func- 
tional networks connecting correlated human brain sites 
[5]. Analysis of the resulting networks in various tasks 
(e.g., move a finger) has shown interesting features in the 
brain, such as scale- free structure, a high clustering coef- 
ficient, and a small characteristic path length. 

Complex network are relevant not only from an aca- 
demical point of view but also from an applied perspec- 
tive. Models based on complex networks for the spread of 
diseases have identified mitigation strategies for epidemic 
spread. In [9] the spread of an infection was analyzed for 
different population structures, ranging from ordered lat- 
tices to random graphs, and it was shown that for the more 
ordered structures, there was a fluctuating endemic state 
of low infection. In [10] it was shown that outbreaks can be 
contained by a strategy of targeted vaccination combined 
with early detection, thus avoiding mass vaccination of the 
hole population. In communication networks, error toler- 
ance and attack vulnerability are key issues. In [llj it was 
shown that error tolerance is not shared by all networks 
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that contain redundant wiring, but is displayed only by 
those with scale-free topology. 

The structure of a network is a key issue in determining 
its functional properties. In real networks communities, 
or modules, associated with highly interconnected parts, 
have been identified [T^[T^ . A nice example of a com- 
munity structure has been unveil in networks of musical 
tastes [13], having also practical applications for the de- 
velopment of commercial music recommendation engines. 

Well-characterized modules have been identified in syn- 
thetic gene networks [12], where positive feedback and 
noise play important roles for the repression and the acti- 
vation of gene expression. An excitable module containing 
positive and negative feedback loops has been identified 
as a key mechanism inducing transient cellular differenti- 
ation [TB]. A method for classifying nodes into universal 
roles according to their intra- and inter-module connec- 
tions has been recently proposed [T7] and applied to the 
study metabolic networks. 

The synchronizability of a network, or is propensity 
for synchronization, is another key issue in determining 
the network functional properties [18] . Why some topolo- 
gies are easier to synchronize than others, is still an is- 
sue not fully understood. Heterogeneity in the connec- 
tion strengths tends to enhance synchronization |19j . A 
weighting procedure based upon the global structure of 
network pathways has been shown to improve synchro- 
nizability |20| . Dynamical adaptation, where the coupling 
strengths develop according to the local synchronization 
of the node and its neighbors, resulting in weighted cou- 
pling strengths that are correlated with the topology, also 
enhances synchronizability [21j . However, heterogeneity in 
the connectivity distribution can have the opposite effect: 
networks with homogeneous connectivity have been found 
to have larger propensity for synchronization than the het- 
erogeneous ones [22] . 

The speed of transmission of information among the 
network also affects the synchronizability. Instantaneous 
interactions have been studied a lot in spite of the fact 
in many situations they are not realistic, because the in- 
formation propagates with a finite speed. A more realis- 
tic scenario considers that the links have associated delay 
times, and that the delays are the same for all the links. 
It has been shown that the presence of such uniform de- 
lays in the communications among the nodes can result in 
enhanced synchronizability. In [23| it was shown that, in 
a network of chaotic maps, it was possible to synchronize 
the delayed network where the undelayed network, with 
instantaneous links, did not. 

An even more realistic approach for complex and disor- 
dered systems is to consider heterogeneously distributed 
delays. In previous work [Ml[^ we studied networks of 
coupled maps with links that have heterogenous delays, 
and investigated the relation of the network topology with 
its ability to synchronize. We found that (i) the synchro- 
nizability was enhanced by random delays as compared to 
networks with uniform delays, (ii) the network synchro- 
nized in a steady state in the presence of random delays 
(in contrast, with uniform delays the synchronization is 



in time a dependent state, [23] ) and (iii) the synchroniz- 
ability depends mainly on the mean connectivity and is 
rather independent of the topology. 

The aim of this paper is to further analyze how the 
synchronizability depends on the connectivity and on the 
topology, when there are delays in the links among the 
nodes. We consider both regular and random network topolo- 
gies, covering the cases of homogeneous and heterogeneous 
distribution of the links. We also study the effects of cen- 
trality and locality on the dynamics of the array. We char- 
acterize the synchronizability in terms of two indicators, 
one that tends to zero when the networks synchronizes, re- 
gardless if the synchronization is in a time dependent or in 
a steady state, and the other that tends to zero only when 
the synchronization is in a steady state. Using these in- 
dicators we also analyze the impact of self-feedback links. 
This paper is organized as follows. Section [5] presents the 
model, the different distributions of delays, and several 
topologies used. Section [3] presents the results of the sim- 
ulations, and, finally. Sec. H] presents a summary and the 
conclusions. 



2 Model 

We consider logistic maps coupled as: 

N 

x^it + 1) = (1 - e)f[x,{t)] + ^Y1 ^^^n^A* ' ^'^^■)]' (1) 

where t is a discrete time index, i is a discrete spatial 
index (i = 1 . . . A^), f{x) — ax{l — x) is the logistic map, 
e is the coupling strength and > is the delay time in 
the interaction between the ith and jth nodes (the delay 
times Tij and Tji need not be equal). The matrix rj = {r]ij) 
defines the connectivity of the network: rjij = rjji = 1 
if there is a link between the iih and jth nodes, and zero 
otherwise. The sum in Eq. (|T]) runs over the hi nodes which 
are coupled to the ith node, 6,; = rjij. The normalized 
pre-factor \/hi means that each map receives the same 
total input from its neighbors. 

A particularly simple solution of Eq. ^ is such that 
all the maps of the network are in a fixed point of the 
uncoupled map, i.e., 

Xi{t)=xi3 V i, (2) 

with xo = f{xo). While this solution exists for all delay 
distributions (rij) and for all coupling topologies (ijij), 
the statistical linear stability analysis performed in |24j 
showed that this solution is unstable unless the distribu- 
tion of delays is wide enough. In the other limiting case of 
all-equal delays. Tij ~ t y i and j, it was shown in Ref. [23j 
that the network synchronizes isochronally, in a spatially 
homogeneous time-dependent state: 

x,{t)=x{t) V i, (3) 

with x{t + !) = (!- e)f[x{t)] + ef[x{t - t)]. 
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2.1 Distribution of delays 

We consider delays distributed as: ~ tq + near(c^), 
where c is a parameter that allows varying the width of 
the delay distribution; ^ is Gaussian distributed with zero 
mean and standard deviation one; near denotes the near- 
est integer. Depending on tq and c the distribution of de- 
lays is truncated to avoid negative delays. 

The synchonizability of the network depends on both, 
the mean delay and the width of the delay distribution 
[24] . The larger the value of tq, the larger the dispersion 
of the delays has to be, for the network to synchronize in 
the steady state (see Fig. 3 of [51]). A similar observation 
was recently reported in [32| . for an integro-differential 
equation describing the collective dynamics of a neural 
network with distributed signal delays. An interesting in- 
terpretation of both observations is provided by the work 
of Morgado et al. [35]: a certain degree of randomness 
in a network (due to random delays, random connectiv- 
ity, or even random initial conditions) results in additive 
and/or multiplicative noise terms in an "effective" single- 
node equation of motion for Xi(t + I), expressed in terms 
of Xi{t) and nonlinear feedback memory terms, that for 
the Logistic map are of polynomial type in Xi (t — n) with 
1 < n < t. In |33| it was found that the right memory 
profile can create optimal conditions for synchronization, 
in other words, an optimal memory range enhances a net- 
work propensity for synchronization. 

Since our aim is to study the roles of the topology 
and of the connectivity, we keep fixed the distribution of 
delays, defined by the parameters tq and c. However, we 
compare the results obtained for distributed delays (c ^ 0) 
with those obtained for all-equal delays (c = 0) and with 
those obtained for no delays (tq = 0, c = 0). 

2.2 Network connectivity and topology 

The connectivity of a network is measured in terms of the 
average number of links per node, 

1 ^ 

w = ]^E^^- (4) 

i=l 

We consider the five topologies, three of them are reg- 
ular networks where the links are distributed determinis- 
tically among the nodes, while the other two are heteroge- 
nous networks where the links are distributed stochasti- 
cally, with given rules. The networks are: 

(i) a nearest-neighbor network (referred to as NN net- 
work) with periodic boundary conditions where each node 
i is linked to its neighboring nodes i±l,i±2,...,i±i4r, 
with K an integer. The number of neighbors is the same 
for all the nodes, and (b) = 2K. 

(ii) the NN network with the addition of a central node 
connected to all other nodes (referred to as ST network). 
In this case, TV — 1 nodes have 2K + 1 links and one node 
has N-l links. Thus, (5) = [(2A'+l)(A^-l)+(A^-l)]/iV = 
{2K + 2){N - \)/N. 



(ill) a star-type network where there are K central 
nodes that are connected to all other nodes and N — K 
peripheric nodes that are connected only to the central 
ones (referred to as KA). In this case, K nodes have 
— 1 links and N — K nodes have K links. Thus, (b) = 
[K{N -1) + {N- K)K]/N. This network has a pure star- 
type structure and is centrally organized^ in contrast to 
NN, that is locally organized. 

(iv) a small-world (SW) network constructed according 
to the Newman and Watts algorithm [55] . 

(v) a scale-free (SF) network constructed according to 
the Barabasi and Albert method [?7] . 

Recently, the influence of feedback loops on the syn- 
chronization of complex systems has received attention. 
In networks of globally coupled units (self-sustained os- 
cillators or maps) delayed feedback in the mean field can, 
depending on the feedback strength and on the delay time, 
enhance or suppress synchronization [28| . In a system com- 
posed by two phase oscillators with instantaneous mutual 
coupling, a delayed feedback loop in each oscillator en- 
hances synchronizability if the coupling strength is not 
too strong, while degrades synchronizability if the cou- 
pling exceeds a threshold value [29] . In small networks of 
mutually delayed-coupled oscillators, which typically do 
not exhibit stable isochronal synchronization, by includ- 
ing delayed feedback loops to the nodes, the oscillators 
become isochronally synchronized [30]l31j. To analyze the 
relevance of feedback loops we consider two situations: the 
diagonal elements of the coupling matrix, rja, are all set 
equal to 1, or are all set equal to 0. In this way, in each 
node a feedback loop is included, or is forbidden. The de- 
lays of the feedback loops, t^,;, have the same distribution 
as the delays of the mutual interactions, r^ . 

2.3 Synchronization indicators 

To characterize the degree of synchronization and to dis- 
tinguish between steady state and time dependent syn- 
chronization, we use the following indicators, 

i 
i 

where (.)« denotes a space average over the nodes of the 
network, (.)t denotes a time average, and xq is the fixed 
point of the uncoupled logistic map, xq = f{xo)- cr^ = 
if and only if Xi = xj \/ i, j, while a'^ = if and only if 
Xj = .To V i. Thus, cr'^ allows to distinguish synchroniza- 
tion in the steady state from synchronization in an time 
dependent state. In the former case, both and a'^ are 
zero, in the latter case, only cr^ = 0. 

3 Results 

In this section we present the results of the simulations. 
We consider a network of A^ = 200 logistic maps with 
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a ~ A interacting with different topologies, as described 
above. The parameters of the Gaussian delay distribution 
are tq = 5 and c = 2. Similar results are found for other 
values of tq and c, with c large enough [24] . To asses the 
role of distributed delays, we compare with two limits: (i) 
instantaneous coupling, = V «, j, and (ii) homoge- 
neous delays, = tq V z, j (c = 0). 

The simulations start from a random initial configu- 
ration, with Xi{0) randomly distributed in [0,1], and the 
maps evolve initially without coupling, during a time in- 
terval < t < max(Tij), because the integration of de- 
layed equations requires the knowledge of the past state of 
the system over a time interval given by the maximum de- 
lay. After that, the coupling is turned on. We neglect tran- 
sient effects disregarding a few thousands of iterations. 

Networks of coupled elements usually show multista- 
bility: different initial conditions lead to different final 
states, and delayed coupling tends to increase the number 
of coexisting states . Multistability is also enhanced 

when the local dynamics of the uncoupled maps presents 
two or more competing attractors [36] , which does not oc- 
cur for the logistic map with a = 4. When the coupling 
is weak and the network is not synchronized, we observe 
a large diversity of dynamical clustered states. However, 
in this "weak coupling regime" the global synchronization 
indicators, cr^ and cr'^, depend mainly on e and only in a 
small extend on the initial conditions, the distribution of 
delays, the network topology and the connectivity (with 
the exception of all-even delays, for them there is a syn- 
chronization "island" that will be discussed below). In the 
following, the plots of cr^ and cr'^ are done by averaging 
over 10 different states generated from randomly chosen 
initial conditions. The only exception is Fig. 6, that is 
done with just one initial condition (that is the same for 
all e). 

In the first section we characterize the network propen- 
sity for synchronization in terms of the indicator cr^ ; in the 
second section, in terms of a'^. In these two sections the 
network does not contain self-feedback links [rja = Q V z). 
In the last section we asses the role of feedback links, by 
studying the same networks (with the same delay distri- 
butions and topologies) but with a feedback link in each 
node (jiii = 1 V i). 



3.1 Characterization in terms of cr^ 

Figures [1] and [2] display color-coded plots of the synchro- 
nization indicator as a function of the coupling strength, 
e, and the average number of neighbors, (6), for the dif- 
ferent topologies and delay distributions. 

Figure[l]displays results for the regular topologies (NN: 
top row, ST: central row, KA: bottom row) and the three 
distributions of delays considered (no-delays: left column, 
homogeneous delays; central column, and distributed de- 
lays: right column). Figure [5] displays results for the het- 
erogeneous networks (SW: top row, SF: bottom row) and 
the same delay distributions. 

We observe a general trend towards synchronization 
when increasing the coupling strength and the average 




0.12 



0.1 



0.08 



0.06 



1 0.5 1o.04 



0.02 



10 0.5 1 



Fig. 1. (Color online) Plot of in the parameter space (cou- 
pling strength, average number of neighbors) for three net- 
works with regular topologies: nearest-neighbors (NN), top 
row, nearest-neighbors with central node (ST), middle row, 
and bottom row K-to-all (KA). The delays are: zero (left col- 
umn), homogeneous {jij — 5 W i and j, central column), and 
heterogeneous { tq = 5, c — 2, right column). 
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Fig. 2. (Color online) Plot of the order parameter a^, Eq[5] 
in the parameter space (coupling strength, average number of 
neighbors) for two heterogeneous networks of 200 nodes: SW 
(small- world, top row), SF (scale- free, bottom row). The de- 
lays are as in Fig.[l] no-delays (left column), homogeneous de- 
lays (center column), and heterogeneous delays (right column). 
Other parameters as Fig. 1. 
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number of links. However we note some important differ- 
ences depending on the topology and on the delays. 

(i) Heterogeneous delays: if (6) is large enough, the syn- 
chronizability does not dependent on the topology, note 
the remarkable similarity of the five panels in the right 
columns of Figs. [T] and [2l for (6) > 10. 

(ii) No-delays and homogeneous delays (left and central 
columns in Figs. [T]and[21): the topology makes a signifi- 
cant difference in regular networks (Fig.[T]), and is less im- 
portant in heterogenous networks (Fig.[2|), although some 
differences can be observed for small (6) ((6) < 10). 

(iii) Regular networks: in Fig. [U comparing the pan- 
els in the right and central columns we notice that the 
network KA (bottom row) is the one with better propen- 
sity to synchronization, while the network NN (top row) is 
the one exhibiting poorer synchronizability. We also notice 
that ST (middle row) has better synchronizability than 
NN when (&) is small, but there are no significant differ- 
ences between them for large enough (b). 

We note that for larger values of (6) (not shown in Fig. 
[1] because we focus on the weak connectivity region) the 
NN network eventually synchronizes. Synchronization oc- 
curs for e above a certain value, e*, that decreases with in- 
creasing (6) , in good agreement with the results of Ref . [37] 
(see in particular Fig. 4(b) of [Sj, where the variable a 
plays the role of (&) here; increasing (&) resulting in a tran- 
sition from local to global couphng). 

(iv) Heterogeneous networks: in Fig. [2] we do not ob- 
serve a significant difference: the panels in the top and 
bottom rows are similar, at least for (5) > 5. 

(v) With homogeneous delays, middle column in Figs. 
[T] and [21 there is an "island of synchronization" in the 
region e ~ 0.15 to 0.19. In spite of the fact that the cou- 
pling is weak, the network tends to synchronize, regardless 
of the topology. This is due to the fact that the delay times 
are even; for odd delays this synchronization window does 
not exist (see, e.g.. Fig. 1 in [25 and Fig. 2 in [36|)- In 
the middle column of Fig. [1] and [21 we observe that the 
"synchronization island" is robust and occurs for all the 
topologies. 



3.2 Characterization in terms of cr'^ 

In order to investigate in more detail the two synchroniza- 
tion regimes (time-dependent and fixed-point synchroniza- 
tion) now we consider the synchronization indicator cr'^, 
Eq. ([S]), that tends to zero only when the network syn- 
chronizes in the fixed point. Figures [3] and [4] display color- 
coded plots of cr'^ as a function of the couphng strength, 
e, and the average number of neighbors, (6), for the same 
delay distributions and network topologies as Figs. 2 and 
3 respectively. 

In the case of heterogenous delays, again we notice a re- 
markable similarity in the five panels in the right columns 
of Figs. [3] and d] for (6) large enough, confirming that the 
synchronizability of the network is largely independent of 
its topology. However, we note that for the topology KA 
(bottom panel, right column in Fig. [3]) the network losses 
synchrony if the coupling is too strong (e « 1). 
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Fig. 3. (Color online) Plot of the order parameter a'^ , Eq. 
in the parameter space (coupling strength, average number of 
neighbors) for three networks with regular topologies: nearest- 
neighbors (NN), top row, nearest-neighbors with central node 
(ST), middle row, and nearest-neighbors with two central 
nodes (KA). Parameters are as in Fig. 1. 



Regular networks: comparing the panels in the right 
and central columns of Fig. [31 we can notice that the net- 
work KA (bottom row) is the one with poorer synchroniz- 
ability in the steady-state. 

Heterogenous networks: in Fig. [4] is observed that the 
SW and SF topologies have very similar propensity for 
synchronization in the steady-state. 

It can be noticed that, as the average number of links 
per node and the coupling strength increase, the networks 
with heterogenous delays tend to synchronize in the steady 
state, while the networks with no-delays or with homoge- 
nous delays do not [in Figs. [3l and [H compare the left 
and central columns with the right column]. This ten- 
dency is particularly clear in heterogeneous networks. Fig. 
[31 where in the left and in the central columns cr'^ is large 
for (6) > 15 and e > 0.8 (right, bottom corner of the pan- 
els) , while in the panels of the right column, a'^ is zero or 
very small in that region. 

The "synchronization island" discussed in the previous 
section is not observed in Figs. [3] and [4] because in this "is- 
land" the network synchronizes in a time-dependent state 
that has a'^ different from zero. 



3.3 Connectivity threshold 

In the presence of heterogeneous delays, it can be expected 
that the average number of neighbors has to be above a 
certain value for the distribution of delays to play an effec- 
tive role. This is indeed shown in Fig.[5l where we plot the 
synchronization indicators and a'^ in the plane (c/tq, 
{b)) for various network sizes and the SW topology, cr^ 
and a'^ are here plotted on a logarithmic scale to reveal 
the following features of the synchronization transition: 
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Fig. 4. (Color online) Plot of a'^ in the parameter space (cou- 
pling strength, average number of neighbors) for three hetero- 
geneous networks: SW (small- world, top row), RN (random 
network, middle row), SF (scale- free, bottom row); and the 
different delay distributions; instantaneous (left column), ho- 
mogeneous delays (center column), and heterogeneous delays 
(right column). Parameters are as in Fig. 1. 



for connectivity below (b) « 10, the network does not syn- 
chronize, regardless of the network size and of the width of 
the delay distribution. Similar results are found for other 
values of e and tq. For (b) > 10, as c increases and the de- 
lays become more heterogeneously distributed there is a 
sharp transition to isochronous synchronization, reveled in 
the plot of cr^ (left column) , which shows a sharp boundary 
between the dark and light regions, at c/tq slightly below 
0.5. The smooth transition seen in the plot of a'^ (right 
column) reveals that initially the network does not syn- 
chronize in the fixed point but in a time-dependent state, 
and as c increases and the delay distribution enlarges, it 
gradually approaches the fixed point. 



3.4 Influence of self-feedback links 

We assess the effect of self- feedback links by setting rju = 1 
V i, in the previously studied network topologies. 

In a large enough network (as the one studied so far, 
with N = 200 maps), the simulations show that feedback 
links have a small influence on the network synchroniz- 
ability when the connectivity is low (when (6) < 0.5A^), 
and their effect is negligible when {b) is larger. With 2D 
color-coded plots is difficult to appreciate the influence of 
the feedback links; therefore, we plot in Fig. [5] the global 
synchronization indicators vs. the coupling strength, for a 
fixed value of (b), that is as low as possible. 

We consider three regular topologies: (i) a ring with 
two neighbors per node and a self-feedback link in each 
node (black solid line); (ii) the same ring without the self- 
feedback links (dot-dashed line, blue online); and (iii) a 
ST network composed by the ring with two neighbors per 
node (without self-feedback links) and a centrally con- 
nected node (dashed fine, red onhne). It can be noticed 
that with no delays and with uniform delays (left and 
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Fig. 5. (Color online) Plot of (left column) and a'^ (right 
column) in logarithmic scale, in the parameter space (normal- 
ized width of the delay distribution, c/tq, average connectivity, 
(&)) for three system sizes: A'^ = 200 (top row), N = 400 (mid- 
dle row), and N = 800 (bottom row). Parameters are ro = 5 
and e = 1, SW topology. 



central columns in Fig. [6|), the ST network synchronizes 
for e large enough, while the ring, with and without self- 
feedback links, does not. With randomly distributed de- 
lays (right column in Fig. [5]), the three networks do not 
synchronize for any value of e. 

Let us now consider the influence of feedback links in 
a smaller network of = 20 maps. Figure 6 displays the 
global synchronization indicator while Figs. 7-9 present 
some examples of the space and time evolution of the net- 
work. The figures are done by representing in a color scale 
the variable Xi{t), with the space index i on the horizontal 
axis and the time index t on the vertical axis. The left col- 
umn of Figs. 6-9 shows results for networks without feed- 
back loops, and the right column, for the same networks 
with feedback loops. A relevant effect of the feedback can 
be observed on both, the global macroscopic indicator and 
on the microscopic network configuration. There are situ- 
ations in which feedback loops enhance coherence, giving 
rise to spatially more ordered patterns (e.g., in Fig. 7, 
top and bottom row), while in others, on the contrary, 
feedback loops destroy the spatial coherence of the pat- 
tern (e.g., in Fig. 9). Interestingly, it can be observed that 
with random delays. Fig. 9, the synchronization is not al- 
ways on a homogeneous state, but there are also static 
patterns with spatial "antiphase" arragement. The char- 
acterization of these patterns is in progress and will be 
reported elsewhere. 
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Fig. 6. (Color online) Plot of (top row) and a ^ (bottom 
row) vs. the coupling strength for a NN network with two 
neighbors per node and a feedback link (solid line), for the same 
ring without feedback links (dot-dashed line, blue online) and 
for the ST network composed by the ring plus a central node 
(dashed line, red online). The delays and other parameters are 
as in the previous figures: no-delays (left column), homoge- 
neous delays (central column) and heterogeneous delays (right 
column) . 
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Fig. 8. (Color online) Synchronization patterns with instan- 
taneous interactions. The network topologies are: a nearest- 
neighbors ring of A'^ = 20 nodes (top row); a A = 20 nearest- 
neighbors ring with the addition of a central node (middle row) 
and a star network composed by a central hub connected to 
N — 20 nodes (bottom row). Without feedback loops (left 
column), with feedback loops (right column). The coupling 
strength is e = 1 
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Fig. 7. (Color online) Synchronization indicator for a small 
network of iV = 21 nodes without feedback loops (left column), 
and with feedback loops (right column). The delay distribu- 
tions are: no delays (top row), fixed delays (middle row), and 
distributed delays (bottom row) . The topologies are: a nearest- 
neighbors ring (black), a ring with the addition of a central 
node (red) and a pure star network composed by a central hub 
connected to A = 20 nodes (blue). 
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Fig. 9. (Color online) Synchronization patterns with homoge- 
neous delays (ro = 5, c = 0). The network topologies are as in 
Fig. 7. (a),(b) e = 1; (c)-(f) e = 0.8. Without feedback loops 
(left column), with feedback loops (right column). 



It is also interesting to notice that when the coupHng 
strength is small (roughly speaking, when e < 0.2) the de- 
pendence of the global synchronization indicators and 
cr'^ with e, shown in Figs. [6] and 6, is very similar. This 
suggests that the global dynamics is almost independent 
of the topology, the connectivity, the network size and the 
delay distribution. We refer to this region as the "weak 
coupling" region. It has been recently shown that in this 
region all the nodes exhibit a qualitatively similar sym- 
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Fig. 10. (Color online) Synchronization patterns with hetero- 
geneous delays (ro = 5, c = 2). The network topologies are as 
in Figs. 7 and 8. (a),(b) e = 0.9; (c),(d) e = 1; (e),(f) e = 0.8. 
Without self-feedback loops (left column), with self-feedback 
loops (right column). 



bolic dynamics, that, for instantaneous interactions, de- 
pends mainly on the network architecture and only to a 
small extent, on the local dynamics [40] , 



4 Summary and Conclusion 

We studied the synchronizability of a network focusing on 
the roles of the connectivity, the topology, and the delay 
times that are associated with the links. The nodes were 
modeled by chaotic logistic maps, and various topologies 
and delay distributions were considered. For low connec- 
tivity (roughly speaking, when the mean number of links 
per node is (b) < O.IA^), and for weak coupling (e < 0.1), 
the network displays an irregular oscillatory dynamics, 
regardless of the network topology and the delay distri- 
bution. For large enough connectivity and coupling, the 
synchronization is mainly determined by the delay times: 
when the delays are homogeneously distributed the net- 
work shows collective synchronous time-dependent oscil- 
lations; when the delays are sufficiently distributed, the 
network synchronizes in a spatially homogeneous steady- 
state. The propensity towards these two synchronization 
regimes was characterized in terms of two indicators, 
and cr'^ [Eqs. ([5]) and the first one tends to zero when 
the network synchronizes isochronously [xi{t) = Xj{t) V 
i, j], while the second one tends to zero only when the 
network synchronizes in the steady state [xi(t) = xq W i, 
with xq — f{xo) being the fixed point of the uncoupled 
maps] . 

When the coupling strength is weak the dependence 
of the synchronization indicators with e is very similar 
for all topologies and delay distributions considered. This 
suggest that in this region of "weak coupling" the net- 
work topology and the delay distribution play no relevant 



role in the dynamics. However, is important to remark 
that the global synchronization indicators employed, 
and ct'^, have the limitation that they characterize only 
isochronal synchronization and steady-state synchroniza- 
tion respectively. In a complex network where the inter- 
actions among the units are not instantaneous, other syn- 
chronization patterns are also expected, such as states 
where the nodes are synchronized but with lag-times be- 
tween them |38il39| . An interesting indicator to analyze in 
future studies is the one that measures the average dis- 
tance between the present state of a map and the delayed 
state of the maps interacting with it: 



{l/N)J2a/b,)J2v^AMt)-^■lit 



Also, it will be very interesting to analyze synchronization 
patterns using symbolic dynamics |40| and complexity in- 
dicators [4T|[42]. 

With heterogeneous delays we also found that there is 
a connectivity threshold below which the network does not 
synchronize, regardless of the network size. This minimum 
average number of neighbors is also independent of the de- 
lay distribution. Above the minimum connectivity level, as 
the width of the distribution increases, the plot of re- 
veals that there is a sharp transition to isochronous syn- 
chronization, while the plot of a'^ reveals that the net- 
work does not synchronize in the fixed point but in a time- 
dependent state, that gradually approaches the fixed point 
as the delays become more heterogeneous. Simulations were 
preformed for SW topology; the study of how this con- 
nectivity threshold depends on the network topology is in 
progress and will be reported elsewhere. 

We studied the influence of feedback loops in each node 
of the network; these feedback loops having the same de- 
lay distribution as the mutual interactions, and found that 
when the network is large, feedback loops have very little 
impact on the global synchronization indicators. However, 
they affect synchronizability of small networks, enhancing 
or degrading the synchronization, depending on the net- 
work architecture and on the delay distribution. As a fu- 
ture study it will be interesting to analyze the interplay of 
delayed feedback with instantaneous mutual coupling and 
vice-versa, instantaneous feedback with delayed coupling. 
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